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Abstract 

We give an exhaustive, non-perturbative classification of exact travelling-wave solu¬ 
tions of a perturbed sine-Gordon equation (on the real line or on the circle) which is used 
to describe the Josephson effect in the theory of superconductors and other remarkable 
physical phenomena. The perturbation of the equation consists of a constant forcing term 
and a linear dissipative term. On the real line stable solutions with bounded energy den¬ 
sity are either the constant one, or of solitonic (or kink) type, or of array-of-solitons type, 
or of “half-array-of-solitons” type. While the first three have unperturbed analogs, the 
last type is essentially new. We also propose a convergent method of successive approxi¬ 
mations of the (anti)soliton solution. 
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1 Introduction 


The purpose of this work is an exhaustive, non-perturbative analysis of travelling-wave 
solutions of the “perturbed” sine-Gordon equation 

- ^xx + sin v? + aipt + 7 = 0 , x ( 1 ) 

for all constant a > 0,'y G M. This equation has been used to describe with a good 
approximation a number of interesting physical phenomena, notably Josephson effect in 
the theory of superconductors P, which is at the base [4] of a large number of advanced 
developments both in fundamental research (e.g. macroscopic effects of quantum physics, 
quantum computation) and in applications to electronic devices (see e.g. Chapters 3-6 in 
0), or more recently also the propagation of localized magnetohydrodynamic modes in 
plasma physics [22]. The last two terms are respectively a dissipative and a forcing one; 
the sine-Gordon equation (sGe) is obtained by setting them equal to zero. 

In the Josephson effect (for an introduction see e.g. Chp 1 in [1]) (p{x,t) is the phase 
difference of the macroscopic quantum wavefunctions describing the Bose-Einstein con¬ 
densates of Cooper pairs in two superconductors separated by a very thin, narrow and 
long dielectric (a socalled “Josephson junction”). The 7 term is the (external) “bias cur¬ 
rent”, providing energy to the system, whereas the dissipative term ai^t is due to Joule 
effect of the residual current across the junction due to single electrons. 

It important to clarify: a) which solutions of the sGe are deformed into solutions 
of ([T|) with the same qualitative features; b) whether ([T|) admits also new kinds of so¬ 
lutions. Candidate approximations to the former can be obtained within the standard 
perturbative method na Eg na El based on modulations of the unperturbed solutions 
with slowly varying parameters (typically velocity, space/time phases, etc.) and small 
radiation components. In particular, the Ansatz for a deformation of a travelling-wave 
solution ip(^Q-^{x,t) = g{Q-^{x—vt) of the sCe reads 

= g(Q){x-xo{t)-v{t)-t^+-fip(i){x,t) + ...\ (2) 

7 plays the role of perturbation parameter, whereas the slowly varying xo{t),v{t) and the 
perturbative “radiative” corrections 7 (^ 7 )(x,t) + ... have to be computed perturbatively 
in terms of a(pt + 7 . If in particular </?(o)(a;, t) is a (anti)solitonic [or (anti)kink] solution, 
one finds PES] also candidate approximate solutions with constant velocity 

i}(t) = Uoo := ±[1 + (4a/7r7)^]“5 ( 3 ) 

which are characterized by a power balance between the dissipative term atpt and the 
external force term 7 . Clearly the perturbative series ([2|) will converge to an exact solution 
only within its ray of convergence, admitted that the latter is nonzero. Even if this 
were the case and one were able to control the convergence, this perturbative approach 
would certainly fail for large 7 . Numerical resolutions [101 E6] of ([T|) are surely a useful 
alternative, but cannot provide solid, exhaustive answers to the two questions above. 

The purpose of this work is to answer questions a), b) by providing (section [3]) a 
non-perturbative classification of exact travelling-wave solutions of ([I]) on the real line 
or on the circle for all a >0,'j G M, and to propose (section E]) an improved method of 
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successive approximations converging to the (anti)soliton solutions, at least for sufficiently 
small 7 . We stick to solutions of physical interest, namely solutions that are stable and 
have bounded energy density h, and therefore also bounded derivatives; in the sequel 
we shall denote them as the relevant solutions. We have begun this job in Ref. [7]. 
The classification is based on a detailed phase space analysis (initiated in [23]) of the 
solutions of the o.d.e. (ordinary differential equation) with varying parameters derived 
(section | 2 |) by replacing in the p.d.e. (partial differential equation) ([I|) a function of a 
suitable combination of x,t only (travelling-wave Ansatz), and on the comparisons of the 
solutions for different values of the parameters; although this seems a rather natural thing 
to do, we have not found in the literature such a classification. If the velocity is ±1 the 
o.d.e. is of first order and can be solved by quadrature, otherwise it is the second order 
one describing the motion along a line of a particle subject to a “washboard” potential 
and immersed in a linearly viscous fluid, and therefore the problem is essentially reduced 
to studying this simpler mechanical analog. A number of useful monotonicity properties 
(section 12.Ih allow us in particular to identify (Theorem 1 in section [3|) four families of 
relevant solutions: three of them (the arrays of solitons for all values of 7 , the solitons 
and the constants only for 7 < 1 ) are deformations of analogous families of solutions of 
the sGe, whereas the fourth family is without unperturbed analog: as each of its elements 
interpolates between a soliton and an array of solitons (see Fig. El), we have called it a 
“half-array of solitons”. The stability of these solutions has been tested numerically. No 
other relevant solutions exist [7| . The families of perturbed solitons and arrays of solitons 
depend on one free parameter less than the unperturbed ones, as the propagation velocity 
V turns out to be a function of a, 7 , which for the soliton coincides, at lowest order in 7 , 
with ([3l) . 

1.1 Preliminary considerations 

Space or time translations transform any solution into a two-parameter family of solutions; 
one can choose any of them as the family representative element. 

The sGe describes also the dynamics of the continuum limit of a sequence of neighbour¬ 
ing heavy pendula constrained to rotate around the same horizontal x-axis and coupled 
to each other through a torque spring [ 21 ]; ip{x,t) is the deviation angle from the lower 
vertical position at time t of the pendulum having position x. One can model also the 
dissipative term —aipt of ([T]) by immersing the pendula in a linearly viscous fluid, and 
the forcing term 7 by assuming that a uniform, constant torque distribution is applied 
to the pendula. This mechanical analog allows a qualitative comprehension of the main 
features of the solutions, e.g. of their instabilities. The constant solutions of ([T]) are 
ip^{x,t) = — sin“^ 7-|-27r/c and y?“(x,t) = sin“^ 7 -I-(2A:-|-Ijvr. The former are stable, the 
latter unstable, as they yield respectively local minima and maxima of the energy density 
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In the mechanical analog they respectively correspond to configurations with all pendula 
hanging down or standing up. 

Our definition of a soliton solution ip is: is a non-constant stable travelling-wave 

solution with ipx,Pt rapidly going to zero outside some localized region. Then mod. 27r it 
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must be 


lim <^(x,t)=—sin ^7, lim <^(x,t) = —sin ^7+2n7r ( 5 ) 

ir—00 x—H-00 

with n £ Z. As we shall recall below, only n = 1 (soliton or kink) and n = — 1 (anti- 
soliton or antikink) are possible [whereas n = 0 corresponds to the constant ips]- In the 
mentioned mechanical model the (anti)solitonic solution describes a localized twisting of 
the pendula chain by 27 r (anti)clockwise around the axis, moving with constant velocity. 
The above condition yields an energy density h (rapidly) going to two local minima as 

+ 00 

X —> ±00. Although this makes the total Hamiltonian H := J h{x,t)dx divergent, the 

— 00 

time-derivative is finite and non-positive: 

00 

H = — J aipjdx < 0. 

— CXD 

[The negative sign at the rhs shows the dissipative character of the time derivative term 
in (JH)]. The effect of 7 7^ 0 is to make the values of the energy potential density at any 
two minima different; this leaves room for an indefinite compensation of the dissipative 
power loss by a falling down in the total potential energy from one minimum to the lower 
next, and so may account for solutions not being damped to constants as t —> 00. 

Without loss of generality we can assume 7 > 0 . If originally this is not the case, 
one just needs to replace (/? ^ —Lp. If 7 > 1 no solutions p) having finite limits and 
vanishing derivatives for x ±00 can exist, in particular no static solutions. If 7 = 1 
the only static solution ip having for x —> ±00 finite limits and vanishing derivatives is 
p = — 7 r /2 (mod 27 r), which however is manifestly unstable. 


2 Preliminary stages of the analysis 


We specify our travelling wave Ansatz as follows: 

i-.=±c-t p{x,t) = g{i)-'K 

^:=-sign(u)^^ p{x,t) = -g{^) 

^:=sign{v)^ p{x,t)=g{C)- 7 r 

^:=x p{x,t)=g{^)- 7 r 

If u = ±1, replacing the Ansatz in ([T]) one obtains the first order o.d.e., 

ag' = -/-smg. ( 7 ) 

We have already argued in [ 7 ] that if 7 < 1 all its solutions yield [ 7 ] unstable solutions of 
dB), except the static constant one p^{x,t) = — sin ^7. The same argument holds also if 
7 = 1 . If 7 > 1 , by integrating one finds 


if u = ±l, 
if v^>l, 
if 0<u^ < 1, 
if u = 0. 






ds 

7 — sin s ’ 
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the denominator is positive for all s G M, so that the solution g is strictly monotonic and 
linear-periodic, i.e. the sum of a linear and a periodic function, so that 

r/s 

g{i + H) = g{0 + 2 vr, E := ^ ( 8 ) 

Jo u(s) 

with u{g) := g'{g) = (7 —sin^r)/®; this will yield (Theorem [T|) a stable solution (p of 
representing an ‘array of (anti)sohtons’ travelling with velocity ±1 (such velocities are not 
possible in the sine-Gordon case). 


In the rest of the section we assume that u 7 ^ ±1. Replacing in ([T]) we find in all three 
remaining cases the second order o.d.e. 

9 ' + +Ug{g) = ? S M, (9) 


which can be regarded as the 1 -dimensional equation of motion w.r.t. the ‘time’ ^ of a 
particle with unit mass, position g, subject to a ‘washboard’ potential energy’ U{g) and 
a viscous force with viscosity coefficient given by 


U{ 9 ) ■= -(cosfif-Fyg), 


a 

\/b“^ - 1| 


( 10 ) 


Note that in equation Q a, u appear only through their combination ( 0111 ) 9 . and that in 
the range |u| G [0,1[ (resp. |u| g]1,oo[) fj.{\v\) is strictly increasing (resp. decreasing), and 
therefore invertible. In Fig. [1] U(g) is plotted for four different values of 7 ; it admits local 
minima (resp. maxima) only if 0 < 7 < 1 , in the points 


9™ := sin ^'y+2kTT, 


(resp. gjf := —sin ^ 7 -|-( 2 A:-|-l) 7 r). 


As 7 —> 1 the points g]J,g^ approach each other, and for 7 = 1 s'™ = g^ = {2k + l/ 2 ) 7 r 
are inflections points. For 7 > 1 no minima, maxima or inflections exist, and Ug < 0 
everywhere. The “total energy of the particle” e := g''^/2 + U{g) is a non-increasing 
function of as e' = —gg'‘^. 

An exhaustive classification of the solutions of equation Q for all values of //, 7 has 
been performed long ago in several works, starting from (see m or [ 2 ] for compre¬ 

hensive presentations). The equation is equivalent to the autonomous first order system 


u' = —fj,u — sing -|- 7 , 

9' = u. 


( 11 ) 


Since the rhs’s are functions of g, u with bounded continuous derivatives, by the Peano- 
Picard theorem on the extension of the integrals all solutions are defined (global existence) 
on all —00 < C < 00 , and the paths [i.e. the trajectories in the phase space {g,u)] do 
not intersect (uniqueness). Each is uniquely identified by any of its points {go,uo). As 
known, the paths may have finite endpoints (limits as ^ ^ ± 00 ) only at singular points, 
i.e. points where the rhs 'sm vanish. These exist only for 7 < 1 and are 

saddles = {g¥,0), nodes, foci or centers = (9^^^)^ 7 < 1 

( 12 ) 

saddle-nodes Ck = {{‘2k + l/ 2 ) 7 r, 0 ) 7 = 1 
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Figure 1 : The potential energy U = 6 —(cos^f+y^f) for 7 = 0 , 7 = . 5 , 7 = 1 , 7 = 1 . 5 . 


Their classification is recalled in Appendix 14.11 Finally, the solutions are continuous 
functions of the parameters //,7 and of (gojUo) (away from singular points), uniformly in 
every compact subset. 

We are going to see that the latter dependences are also monotonic. To analyze them 
and the asymptotic behaviour of the paths near the endpoints it is useful to eliminate 
the ‘time’ ^ and adopt as an independent variable the ‘position’ g, as in the unperturbed 
case. The path of any solution g{^) of ([ 9 ]) is cut into pieces by the axis u = 0 . Let 
X =]^_,^+[C M be the ’time’ interval corresponding to a piece, 

e := sign(u(0), ^^X 

be its sign and let G =]g-, g+[:= g{X). In X the function can be inverted to give 
a function ^ : g G G ^ ^{g) G X. So one can express the ‘velocity’ u and the ‘kinetic 
energy’ z := tt ^/2 of the ‘particle’ as functions of its ‘position’ g. By derivation we 
find that g”{C) = Ug{g{^))g 'and the second order problem ([ 9 ]) with initial condition 
(fl'(Co)) ^^(?o)) = (sO) in X is equivalent to two first order problems: the first is 

Zg{g) + egy/z{g)+smg--f = uug{g) + gu{g) + smg--/ = 0 , u{go)=uo ( 13 ) 

(note that this is invariant under the replacement 5 —> 5 + 27 r), which has to be solved 
first, and yields a solution u = u{g; go,uo; continuous in all arguments (away from 
singular points); the second is 

g'iO = u{g{^)), g{^o) = go, 

is integrated out by quadrature 

Co 90 go 


( 14 ) 

( 15 ) 
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and implicitly yields a solution g = g{^] gQ,uo; n,'y) in X. If X is not the whole M, the 
final step is the patching of solutions in adjacent intervals X. 

Choosing in (jlSp g as g± one obtains ^±. If z{g) vanishes as ■=\g±—g\^ with a > 2 
as g ] g+ oi g I g-, then = oo or = —oo. The behaviour of u{g), z{g) near g± can be 
determined immediately solving (jl.Sp at leading order in a left (resp. right) neighbourhood 
of g+ (resp. g-). In particular, if 7 < 1 and g± = g^ (a maximum point of U) then the 
equation obtained by replacing the power law Ansatz u{g) = in (jl 3 p is 

solved by 


“(s') ~ ( 5 '+- 5 ')“+e as 5 T 5 +, 

“( 5 ) ~ {g-9-)u-e as g i g-, 


( 16 ) 


where for e, e' G {+, —} Ug/g is defined by 


:= 


:= ^ +4y/’l-72^ . 


Formula (I 16 h gives the leading behaviour of the four separatrices having an end on A^- 


Problem (I 13 p is also equivalent to the Volterra-type integral equation 

9 

z{g) = zo + U{go)-U{g)- e j ds^i^/ 2 z{s) ( 17 ) 

90 

where zq := ng/ 2 . When g = 0 (no dissipation) this gives the solutions explicitly and 
amounts to the conservation of the total energy of the ‘particle’ e{g) = z{g) + U{g). 


2.1 Monotonicity properties 

In agreement with the physical intuition, the solutions of (I13|l and the extremes of 
G depend on the parameters g,zo,^ monotonically. In the Appendix we prove the 
following propositions, which include and extend the results of m^- 

Proposition 1 As functions of zq: z = tt^/2 is strictly increasing; g^ is increasing and 
g- decreasing (strictly as long as they have not reached the values ±00). 

Proposition 2 As a function of both g,—e'y the solution u{g; go,uo] g,'y) is strictly de¬ 
creasing (resp. strictly increasing) for g G\go,g+[ (resp. g G]g-.,go[). Correspondingly, 
the solution g{f,; go,uo; g,'y) is strictly decreasing as a function of both e/r, —7, and so is 
either extreme g± (strictly as long as it has not reached values ±00). 

Remark. In general g± will be discontinuous functions of g,zo,'y at g± = gjf. 

Whenever the domain G of the solution z{g) contains a whole interval ]g,g-\-2TT[ we define 

S+27r 

Hz,g) ■= J dsy/2z{s) ( 18 ) 
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Given any go G G, let gk ■= go+2TTk, K := {kGl^ \ gkGG} and Ik ■= I{z, gk) if k, k-\-\ G K. 
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Proposition 3 If e = — the sequences {z{gk)},{Ik} are strictly increasing and diverging 
as k ^ oo, with K bounded from below. If e = + the sequences {z{gk)},{Ik} are: either 
constant, with K = Z; or strictly increasing and converging as k ^ oo, with K bounded 
from below; or strictly decreasing, diverging ask —oo, and either converging as k ^ oo, 
or with K upper bounded. Moreover, 

z{gk+i) - z{gk) = 27r7 - egh- (19) 


3 The relevant solutions 


li^ = a = g = t) (sGe) the total energy e of the ‘particle’ is conserved and its value 
(together with the free parameter v) parametrizes different kinds of solutions of Q [3l[6]. 
Plotting U{g) (Fig. [2] left) we get an immediate qualitative understanding of them. They 
all have bounded 2 ; and therefore g'. This implies that also the corresponding ipx,(pt,h 
are bounded functions of x,t. 

If 7 , a, /i 0 (perturbed sine-Gordon) there are [3 solutions g{f,) with g' (and therefore 
also (fx, h ) diverging at space and time infinitjo; as said, they are to be discarded. In 
Ref. [7] we have analyzed all the possibilities for 7<1 and shown (Prop. 1) that relevant 
(in the sense of the introduction) solutions ip, if they exist, can be only of four types, all 
with < 1 and e := sign(g'') > 0; out of them three are deformations of travelling-wave 
solutions of the sGe. Actually a simple inspection shows that the arguments and the 
conclusions given there continue to hold for 7 > I; the analysis could be simplified with 
the help of the monotonicity properties of section [2Tl Step by step, we are now going to 
see (Theorem [T]) that all four types actually exist. 

We have already given in the introduction the constant solntions (p^,(p'^. They 
correspond to g{f,) = g^,g]f respectively, implying that e takes the value of a maximum 
or minimum of U{g) (e = ±I if 7 = 0), respectively. 

If 7 = a = /i = 0 solutions corresponding to e = ei g] — I, I[ are unstable [HIEIE]; their 
deformations for 7 , 07 ^ 0 are not relevant as well [7]. 

The (anti)soliton solutions are obtained from non-constant solutions g{ff) having finite 
limits at both ^ ^ ± 00 . The corresponding paths are heteroclinic orbits ending at two 
neighbouring saddle points, e.g. Ao,Ai: the ‘particle’, confined in the interval g^ <g< 
g^, starts at ‘time’ = —00 from one extreme and reaches the other one at ^ = 00 (the 
extremes being maximum points of U{g)). The corresponding ’kinetic energy’ z{g) will 
be defined in the same interval and fulfill the boundary conditions 


lim z{q) = 0, 
gl9^ 


lim z(q) = 0. 


( 20 ) 


If 'y = a = g = 0 these are characterized by e = e = I, so that g{^) —> Tvr as ^ > ± 00 . 

Explicitly, these two solutions can be obtained from (fTTll . (fT^ inverting f,{g). Replacing 

^For instance, by Prop. [3] if e = — and z{g) is defined at least in an interval of length 27r then g_|_ = 00 , z{g) 
diverges as g ^ 00 , g'{f), Px, Pt diverge as ^ ^ — 00 . 
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the result in Q, mod. 27 r they translate into unstable solutions of the sGe if > 1 (all 
pendula stand upwards outside a small region), and the celebrated families 

— vt 
\/l— 

of stable [HEIE] solutions if <1\ all pendula of the chain hang downwards outside 
a small region travelling with velocity u; within that region they twist n = ±1 times, 
i.e. once clockwise or anti-clockwise, around the x-axis. is the one-parameter 

family of soliton solutions, (p'^Q^{x,t]v) the one-parameter family of antisoliton solutions, 
the parameter being their velocity u, which can take any value in ] —1,1[. 

If 7 > 0 an heteroclinic orbit exists only if 7 < 1 , € = + and ^ is adjusted to a special 
value /i(7)>0. In fact, consider equation (fT 3 ]l with varying ^ > 0 and impose (pUji i . For 
/r = 0 the total energy e is conserved: the corresponding solution z{g) is therefore defined 
in gi := gQ + 2 T: = g^ and z{gi) = 2 irj > 0 , whereas for sufficiently large g [ 23 ], g > g, 
the corresponding z{g) is not defined in gi (as g- < gi). By continuity in g and Prop. 
[21 there exists a unique //(y) >0 such that the corresponding solution z{g)) is defined in 
G =]g^,g^[-, i-e. fulfills also (1201) 9 [see Fig. [ 2 | right, where also e{g) is plotted]. By (IT^ 
relation (I 20 h implies e = -|- and 

A-^(2,5^) = 27r7, (22) 

which is the energy balance condition between the energy dissipated by the viscous 
force and the potential energy gap after a 2 t: displacement of the ‘particle’. Replacing 
z{g) in (11511 one finds the mentioned g{^) (this will have g' > ^ everywhere); replacing the 
latter in ([6]) one finds deformed (anti)soliton solutions, as described in Theorem [1] 
By inversion of (110^ 9 the velocity v will no more be a free parameter, but a function (in 
the range ] —1,1[) of 7,0 and of the helicity ± of the (anti)soliton solution 0 ^. 

If 7 = 1 also a path connecting Cfc, Ck+i can be found but this will be of no use 
because it yields a manifestly unstable (p. 

The “array of (anti)soliton” solutions are obtained from linear-periodic solutions g, 
in the sense ([8|): the ‘particle’ travels towards the right from g- = —00 to (7+ = 00 (or 
towards the left from g^ = 00 to g- = —00) and its ‘kinetic energy’ z{g) is 27 r-periodic, 
in particular takes the same value zm at all points g^, 

Kgk) = ^M VA: e Z. ( 23 ) 

If7 = Q: = ^ = 0 the corresponding solutions 5(0) (=*=?; g) can be equivalently characterized 
by the constant e = e > 1 , or by zm = ^ — U{g^), or by the period [because of 
([8|)2]. Again, the corresponding solutions of the sGe are [HEIEI unstable if > 1 and 
stable if < 1 (‘most’ pendula down in the pendula chain model). The stable solutions 
0 '^Q^{x, t) = 5(0) (i?) respectively describe two-parameter families of evenly spaced “arrays 
of solitons and antisolitons”, the two parameters being their velocity u, which can take 
any value in ] —1,1[, and one of the variables e,ZMjS(o)- 

If 7 > 0 , solutions fulfilling ( 1231 ) exist only if e = -|- and g is adjusted to a special value 
g{'y,ZM)> 0 . In fact, consider problem (IT^ with varying g, zq = zm>Q and set 

(9^ 

90 = < ^ 

[ 2A:7r— ^vr if 7 > 1 , 
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Figure 2: The potential energy U{g) = 6 — {cosg + jg) for 7 = 0 (left) and 7 = .1 (right). 
Correspondingly, the ‘kinetic energies’ and the ‘total energies’: 1) z,e associated to the soliton, 
/i = il{j); 2) z,e associated to an array of solitons, p< fj,{j); 3) z,e associated to a half-array 
of solitons solutions, g,<fi{'j). 


Because of (fTF|l condition (125]) implies again e = For p = 0 the total energy e is 
conserved: the solution z{g) is therefore defined in gi := go+2'K and z{gi) = zm + 27r7 > 
zm, whereas for sufficiently large p [25], p > p, either it is z{gi) < zm or even z{g) is 
not defined in gi (as < gi). By continuity in p and Prop. [21 there exists a unique 
g{l:ZM) £] 0 ,p[ such that the corresponding solution z{g) := z{g; go, zm] g,'y) fulfills the 
condition (125]) [see Fig. [2] right, where also e{g) is plotted]. That this holds not only 
for k = 0, but for all k follows from iteration in successive intervals. Actually, by the 
invariance of (|13p under g —> g+2TT and the uniqueness of the solution this implies that z 
is defined in all M and periodic: 

z{g) = Hg+27r). (24) 

By (119]) this implies 

h^(2,p) = 27r7; (25) 

the Ihs is independent of g (and can be called simply I). Again this equality amounts 
to an energy balance condition: the energy dissipated by the viscous force equals 
the potential energy gap after a 27r displacement of the ‘particle’. For g fixed z, I are 
strictly increasing, continuous functions of zm by Prop. mm whereas fi and H are strictly 
decreasing and continuous respectively by (]25l) and ([ 8 |) 2 . All these functions are therefore 
invertible, and one can adopt any of the four parameters p,H (in the appropriate 

range) as the independent one, beside 7 . For |t;| < 1 one can adopt also juj as the 
independent parameter, as the function g{\v\) defined in ([lOh o is strictly monotonic. 

Replacing z{g) in (|15p one finds the mentioned g{^), and replacing the latter in ([ 6 P 
one finds one-parameter families of evenly spaced “arrays of solitons” and of evenly 
spaced “arrays of antisolitons” , as described in Theorem [H as a parameter one can 
choose Zm, I, j^’j- The (anti) soliton solution is recovered as the zm^O limit. 

As we are going to see in the next theorem, the z{g) are asymptotic solutions [23] at 
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g oo, more precisely exponentially attract all other solutions z(g) of eq. (I13p (with the 
same value of /r) with domain extending to g^ = oo, in particular, all z such that z> z] 
the latter fact holds also if /i = /l and we denote by z the “spurious” periodic solution, i.e. 
the periodic extension of z. The z{g) such that z<z will yield a family of solutions of ([l]) 
without unperturbed analog, the “half-array of (anti)solitons”. 

We have partly proved and are ready to state the following theorem, which includes 
and completes results partly obtained in Ref. [7]. 


Theorem 1 Let 


x—vt 


v{g) := 




< 1 . 


(26) 


Mod. 271, travelling-wave solutions of HP (where 7 > 0 and a> 0) having bounded 

derivatives at infinity are only of the following types: 

1. The static, uniform solution ip^{x,t) = — sin“^ 7 , if < 1. 

2. The soliton fi^{x,t) = g{f,) — 7T and the antisoliton (p~{x,t) = g{—ff) — 77, 

lim (,p^(x, t) = — sin“^7, lim (,3^(x, t) = — sin“^ 7 ± 27 r, (27) 


travelling respectively rightwards with velocity v = v := v{fi{'y)) and leftwards with 
velocity v = —v, only i /7 < 1 and a > 0. The function fulfills the bounds |R7P , 
is continuous and strictly increasing in [ 0 , 1 ], with /i( 0 ) = 0 ; a good numerical 
determination is in \24^ , whereby ~ 1,193. It can be determined with arbitrary 
accuracy by the method described in Thm. which gives = 7 r 7/4 -|- 0 ( 7 ^). 

3. The “array of solitons” {x,t; g) = g{f^]g) — 7T and the “array of antisolitons” 
tp~{x,t;g) = g{—(^] g) — TT, only if a > 0, and for any g e] 0 , 00 ] i /7 > 1 , for any 
g G] 0 ,/i( 7 )[ i /7 < 1, where g fulfills 

gPl-K 

r (j c 

g{i + H) = m + 277, E{g,j) = / e] 0, oo[ (28) 

J V2z(s) 

9 

(g is “linear-periodic”), travelling respectively rightwards with velocity v = v{g) 
and leftwards with velocity v = —v{g). //7 > 1, (p^{x,t-,oo) := g{Mr — t) — 7r = 

lim^_,oo g ^ ^ 73 ^* ; the g found in ( 0 ) (these have velocity u = ±1 respec¬ 

tively). 

4 . The “half-array of solitons” (:p~^{x,t) = g{f,]g,3) — T^ and the “half-array of anti¬ 
solitons” (p~{x,t) = 5 (—/;/r, 7 ) — TT travelling respectively rightwards with velocity 
V = v{g) and leftwards with velocity v = —v{g), only i /7 < 1, a > 0, and for any 
g G]0,/i(7)[, where g fulfills g<g and 

lim (^^(x, i) = — sin“^ 7 lim [z(fi() —^( 5 ')] = 0“ (29) 

a;—».i|ic>o S'— 


lim b(0-5(0] = lim b'(0 - 5 (01 = 0 (30) 

^—>00 

^The stability has been tested numerically; an analytic study will be done elsewhere. We just note that the 
key property used in the stability proof of [19], g'{f)>0 is fulfilled by the families of solutions 2,3,4. 
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(for an appropriate choice of g in the family of g’s differing only by a (^-translation). 
The last three limits are approached exponentially. 

All g,g,g,g — g are strietly increasing. To parameterize the solutions of classes 3,4 one 
ean adopt as an independent variable alternative to jj. either zm,I, 1^1 or H. 


Rest of the proof: In Ref. [7] we showed that solutions S' of ([9]) yielding solutions (p 
of ([I]) which have bounded energy density h and are not manifestly unstable can only be 
of the above type. We have just shown that the first three families of solutions actually 
exist, and added some details to their properties. Formula (I26p follows from inverting 
(jl 0 p 9 in the branch |u| < 1 (whereas the branch |u| > 1 yields no relevant solutions). 
Let g(q) := g{(() with g := = sign(u)x — |u|t. gri,grjri are periodic. For 7 > 1, 

replacing in Q and letting |u| | 1 we find that g fufills ([7]). This proves the limit 

lim^^oo 9 ^~ 9{'^~t)i after noting that by (fT 0 |l 9 /r —> 00 as |t;| t 1 - 

As for the fourth family, consider a z{g) with domain extending lo g+ = 00 . Since 
the two diagrams z{g),z{g) do not intersect, w{g) := z{g) — z{g) is either positive- or 
negative-definite. By (fT^ it fulfills 


Wg = -g, 


y/2{J+w) — 


= —2g 


w 


\/ 2 (z-|-'U?) -\-^/2z 


implying 


d . . 2g g 

— In kc =- , -= <- , 

dg y/2{z+w)+V^ ^,/2{z^ + \w{go)\) 


(we have denoted by z^ the maximum of z): |rc( 5 ')| is strictly decreasing. By integration 
we find for g > go 


w^(5)l < ^(50)16 


C := 


_ If _ 

y/2{z^ -P\w{go)\)’ 


(31) 


namely |rc(( 7 )| ^ 0 exponentially as g ^ 00 , as claimed. 

In the case 7 < 1, /U < g{'y) this applies in particular to the solution z fulfilling the 
initial condition z{gjf) = 0 (for some A: G Z). Since g < g{'j), by Prop. [2] this is defined 
and larger than z, therefore positive in g^^ and then by Prop. [3] will be defined in all 
G =]g^,oo[. Since Ak = {g^,0) is a saddle point, the corresponding path {g{(),u{()) is 
a separatrix and 

giO 9]f, 5(0 oo- (32) 

As z{g^) > 0, the corresponding re = z —z is negative-definite, and we find in the order 

w:=z-z\0, V^-V^\0, —^-^ i 0, (33) 

exponentially as g ^ 00 . The first limit gives (|29l) 9. From m 




ds 

—, ~\~c. 

a/2IM 
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where c, c are integration constants, whence 


Cig) - ^{g) = [ ds 
J go 


V^ 2 i(s) ^/2z{s) 


+ (c-c). 


The integrand is positive and goes exponentially to zero as g' —> oo, therefore the integral 
converges. Choosing the constants so that ds[] = (c—c) we find 


^(g) = i(g) -p(g), 



1 

\/2z(s) 


with p(g) positive and exponentially vanishing. Applying the inverse g(^) of ^(g) to both 
sides we find 

g = g (i(g)-p(g)) = g {l{g))-g'{l)p{g)- 

The second equality is based on Lagrange theorem, where ^ is a suitable point in ]'^( 5 ) — 
P{g)d{g)[- Finally, setting g = g{i) we find 

g{i) = g{i)-g{i)p{g{0), 

where ^ G].^ —/ 5 ( 5 (.^)),^[. The second term at the rhs exponentially vanishes as ^ oo 
[since p{g{C)) does and g' is bounded], proving (|3n]) i . By ([30]) i the second limit in ([33]) 
now implies (|30|l 9. 

□ 

Remark [S].! In Fig. [3] we have plotted a soliton, array of solitons and “half-array 
of (anti)solitons” solutions, determined numerically. The latter has no unperturbed ana¬ 
log. It interpolates between the (anti)soliton solution at one extreme and the array of 
(anti)solitons solution at the other. Therefore it cannot be approximated, nor can it even 
be figured out, by the modulation Ansatz (j 2 |). 

Remark [31.2 Setting X := H\/l — v'^, eq. ([28]) and (| 6 ]) imply 

(p^ {x+X,t) = (p^{x,t) ±27r; (34) 


so makes sense also as a solution of ([T]) on a circle of length L = mX, for any m G N. 
The integer m parameterizes different topological sectors: in the m-th the pendula chain 
twists around the circle m times. 


Remark [3l 3 We emphasize that, in contrast with the unperturbed soliton (and array 
of solitons) solutions, where v was a free parameter of modulus less than 1 , u is predicted as 
a function of 7 , a for the perturbed soliton, as a function of 7 , a and one of the parameters 
for the perturbed (half-)array of solitons. 

We determine the ranges of the various parameters. Clearly, as zm —> 00 z and I 
diverge, whereas fi, H, v go to zero. 

We now consider the limit zm ^0. If 7 > 1, as zm —> 0 one finds the following leading 
parts and limits 


7-1 


/Zm 


00 , 


00 , 


7-1 


7-1 


^22m«^+(7-1)^ 


0 ; 


(35) 
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Figure 3: The soliton, the array of solitons and the half-array of solitons solutions 


as zm spans ]0, oo[, the range of any of I, p,, H is ]0, oo[ and that of v is ]0,1[. In fact, the 
Taylor formula of second order for z{g) around can be written without loss of generality 
in the form 

zig;zM-,^) = ZM + ZMC,i{zM-,i){g - gu) + {g- gkfp{g) (36) 

with p{g) bounded; in order that, as zm 0, z keeps nonnegative both in a left and a 
right neighbourhood of gk, Ci{zM','y) has to approach a finite limit. Replacing this Ansatz 
in m we find at lowest order in {g — gu) 

zmCi + PV^zm -|- 1 — 7 = 0. 

As Zm —> 0 this implies (l3^ [by (| 2 ^ . (I28]l 9 and (I 26 II 9 ] . Summarizing, as zm spans ]0, oo[ 
the range of any of I, p, E is ] 0 , oo[ and that of v is ] 0 , 1 [. 

If 7 < 1, by the monotonicity property p{'y, zm) < and by the continuity we find [23] 

lim /i( 7 , zm) = Kl) < 00 . 

ZM^O 

Hence if 7 < 1 the range of p as zm spans [0, oo[ is ]0, p], the range of I is ]2TT'y/g,oo[ the 
range of v is [0,v[. The following bounds for p{'y) have been proved ([23l[9], see [l7] for 
a summary) 

^+ < A( 7 ) < ^ 2 ( 1 - 71 ^)., (37) 
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4 Method of successive approximations 

Eq. (fTTl) can be reformulated as the fixed point equation 

Az = z (38) 


for z{g), where for e > 0 the operator A = A{gQ, zq] g, 7 ) is defined by 

9 


Aw{g) := aj{g-,go,zo-,-f)-Jds(l){g,s,w{s)) 

90 

‘^{g]9o,zo-,'y) ■■= zo + U{go)-U{g) 4>ig,s,C) ■= '/Kg 


(39) 


on the space of nonnegative smooth functions tc on M (the domain of w can be always 
trivially extended to M). According to the method of successive approximations, after a 
reasonable choice of a function 2 ;(o)(fl') as an initial approximation for z{g), better and 
better approximations should be provided by := A"’Z(o)(5') as n ^ 00 . For this to 

make sense, at each step it is necessary that belongs to the domain of A (in the 

present case, it must be nonnegative, otherwise the integrand function is ill-defined) and 
that the sequence converges. With the known standard theorems, this can be guaranteed 
a priori not in the whole domain G of the unknown z, but only in some smaller interval J 
containing go. In general only the iterated application in infinitely many adjacent intervals 
allows to extend a local solution to a global one, what makes the procedure of little use 
for its concrete determination. 

Estimating the length of such a J one finds that it is not less than 27r only for suf¬ 
ficiently large zq. Actually, the determination of the solution in an interval of length 27r 
would be enough for the complete determination both in the case of a periodic solution 
i (which is then extended periodically) and of a separatrix z (in that case G =\g/i,g^[, 
which has exactly length 27r). The periodicity condition (1231) is automatically fulfilled by 
each Z(^n) if we modify the definition of A adjusting the coefficient /x to rc as follows: 


Aw.= A{go,zo;fi{w),'y) w 


g{w) := 27r7 


■ gcrl-27r 

J ds\/2w{s) 

- go 


(40) 


Choosing go = g^^ for simplicity, then ii{z(^n)) will converge to g{'y,zo). If instead we fix 
g, as an independent parameter, one will obtain zo as Z(j{^{go) [23] • For the periodic 
solution a sufficiently large zo amounts to a sufficiently small /x; in [23| the following 
quantitative condition was found: 


hi > ei, 


27rV2 


(41) 


where 

ei := max|z(i)-2;(o)| = lk(i)-^^(o)Iloo, :=min|z(i)|. 

So ? 7 i cannot be too small, in particular cannot vanish, what excludes the cases of the 
periodic solutions i having low energy and of the heteroclinic orbit i. 
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4.1 The soliton solution by the method of successive ap¬ 
proximations 

The standard theorems fail for i because the sup norm has not enough control to guarantee 
non-negativity of the approximations z^n) everywhere in G, as well as the fulfillment of a 
Lipschitz condition by the integrand cj) and the behaviour (jl 6 [) near the extremes of G. In 
this section we adopt a tricky, nonstandard choice of the norm and show (Theoreml^ that 
a single application of the method of successive approximations gives the soliton solution 
{fi,z{g)) in its whole domain G =\g^i,g^[- 

Assume 7 < 1. Choose g^ = g^i, -^o = 0 and \et y ■.= g — go- Then 

uj{y) = y/l- 7 ^ 2 sin^ | + -f{y-s\\iy) = ^\/l-7^2/^ + 0{y^) (42) 

and z fulfills (j38p . where the operator A has taken the form 

y 

Az{y) = z{y)-.= y/l^^2 sin^ |-k 7 ( 2 /-siny)-/ 2 (z) J dy'yj2z{y'), 

0 

27r^ 

where fi{z) := — - (43) 

f dy'^2z{y') 

0 


By (fT6l) z{y) = O(y^), z{2'K—y) = 0^(27r—7/)^j . One easily checks that, more generally, 

if z has such a behaviour near 0, 27r, so has Az. So it would be more natural to look for 
the solution from the very beginning in a functional space whose elements have such a 
behaviour. In (^^([O, 27r]) introduce the norm 


sup 

ye[ 0 , 27 r] 


My) 

p2(y) 


(44) 


where the ‘weight’ p should vanish as y and 2TT — y at 0,27r and will be specified later. 
Clearly 

Ikll > C'lkiloo = C sup \z{y)\ G~^ := sup ^ . (45) 


The subspace 


j/e] 0 , 27 r[ ye] 0 , 27 r[ 

V := {z{y) G C'°°([0, 27r]) | || 2 ;|| < 00 } 


(46) 


is a complete metric space w.r.t. the metric induced by the above norm. In fact, consider 
a Cauchy sequence {zn} C C in the norm || • ||: by it is Cauchy and therefore converges 
to a (uniformly continuous) function z{y) also in the norm || • ||oo; moreover for any e > 0 
there exists f G N such that Vr > f, Vm G N 


sup 

ye[ 0 , 27 r] 


Zrjy) - Zr+m{y) 

p^{y) 


e 
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Letting m ^ oo we find 


sup 

2/e[0,27r] 


Zriy) - z{y) 


p^{y) 


< e, 


showing that z G Vg and that {zn} —>■ z also w.r.t. the topology induced by the above 
norm. 

Let a, 6 E M with b > a > 0. The subset 


Za,b,P ■■= [z{y) G y I < ^ < 62| 


(47) 


is clearly closed w.r.t. the metric induced by the above norm. We shall look for (z, jl) 
within a suitable Za^b,p- First we look for a, b such that (1431) defines an operator A : 
Za,b,p Zafi^p. Up to a factor, we choose p^{y) as the 7 = 0 (i.e. unperturbed) soliton 
solution zo{y), more precisely p{y) := sin |. Then 


2n 


P{y) := J dyp{y) = 2(1 - cos ^) = j dyp{y), 

0 y 

and, since 1 — \/l — w > w/2 we find (setting w = sin^ |) 

P^{y) < P{y) < 2 (^ 1 -cos I) (^ 1 +cos I) = 2 /(y). 
Thus for any 2 ; G Zafi^p we find 


aP{y) < j dy'^2z{y') = J dy' \ {y') < bP{y) 


2-k 


4a = aP(27r) < = f dy^/2z{y^) =< bP{27r) = 46 

M J 


implying the inequalities 77r/26 < jl < 7'7r/2a and 


y 

7^P^(y) <pj dy'yj2z{y') < 'y^p^{y). 


Similarly, 


2-k 


7^/(y) <pj dy'^2z{y') < 7 ^/( 2 /)- 


(48) 


( 49 ) 


Lemma 1 For all y >0 

y‘2. y3 

1 — cos y > 0, y — sin y > 0, - 1 + cos y > 0 ,- y + sin y > 0. 

2 6 


^If per ahsurdum sup \z/p^\ = 00 then the Ihs would certainly exceed e. 
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Proof: The first equality is obvious; the others are obtained by iterated integration over 
[0,y]. Q.E.D. 

As a consequence, for ye [0, vr] 


0 < y-siny < — = 
b b 


1 2 


Sin 


2 J 


y <P^(y)y. 


Collecting the results, on one hand assuming 1 > a/b > 1/2 we find 


z{y) >p^{y) 


21/1—72 — 77 r— 


>y^(y)2 1/1-72 - 77r 


for all y e [0, 27r]; on the other hand, for y e [0, vr] we find 

r3 


z{y) <p (y)2 


-/I-72 + 7 


TT 

T2 


(50) 


(51) 


(52) 


This provides bounds for y G [0, vr]. To find bounds for y G [vr, 27r] set v = {27: —y) and 
note that from (I43p it follows 


z{y) = y/l—72 2sin^ ^ —7(1; —sinu) + 27r7 — 




27r 27r 


J dyV^- J dy'^/2z{y') 


2tt 


= 1/1-72 2sin^ |-7(u-sinu)+/i j dy'^/2z{y'), 


We use to bound the third term at the rhs; as u G [0, vr], to bound the second term 
we can use ([50]) with y replaced by u, but keeping p^(y) = p^(u) at the rhs of the latter. 
Collecting the results we thus find for y G [vr, 2vr] 


y/y)2 y/1-7^ - <z{y)<p^{y)2 + jt: 

Hence <2z< b'^p^, so that z G Za^b,p, if we define 

yl —72 — 77r 


:= 4 

In order that 1/2 < a/b it must be 

1 ^ -/l-72 — 77r 

4 52 ./l —72 -)- ^7]- 

which gives, after some computation. 


b^ := 4 y/l —72 + 77r 


7 < 


1 + 


25vr2 


.187 


(53) 


(54) 


(55) 


We conclude that in this 7-range with the above choice of a, b AZa^b,p 7 Za^b,p, as required. 
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Let us determine the constraints on a, b following from the condition that A be a 
contraction. First, we immediately find 

0| r ^ ^ M 2. ■,‘^\zi{y) - Z2{y)\ . 2( Ml II 

2 \zi{y) - Z2{y)\ = p [y) --< p {y)\\zi - ^2 

p\y) 

Note that for any a > 0, \y/ui — y/v^\ < |rti — U2\/{2a) if ui,U 2 G [a^,oo[. Hence 


\^/ 2 zi{y) - ^/2z2{y)\ = p{y) 



< 


p{y) 

2 a 


Hi - Z2 



^ p(y) ‘^\zi{y) - Z2{y)\ 
2 a p^{y) 


(56) 


1-1 -1 I 

\pi - m = piP2\Pi - P2 1 -^ 


2-k 


dy{yj 2 zi{y) - y/ 2 z^ 


27r 


2tt 


< ^ I dy y/ 2 zi{y) - ^/2z2{y) 


7r7 


j^vp{y) 


Try 

4a^ 


\zi - Z2I 


(57) 


h-zi = j dy' {pi - p2)^/2zl{y') + p2{\/2zi{y') - i/2z2{y')) 
0 


whence 


Hi 


y I - y 

{y)-Z2{y)\ < \PI-P2\ J +^2 J dy' i/2zi{y')-i/2z2{y') 

0 * 0 


- “ ^2||L’(y) + ^IIh - z2\\P{y) 


< 1 + 


Try 

4a? 


ki - z2\\P{y) < 1 + 


Try 


ki - Z2\\p^{y), 


implying 


\h{y) - Z2{y)\\ < ( 1 + -^ ^IHi - Z2ll- 


a / a^ 


Thus, H is a contraction if 
that is, 

7 < 


A ;= (1 + 6/a)Try/a^ < 1, 

,2 4 


< ^ = - 77r], 


-(1 + ^) 


Stt Stt 


(58) 

(59) 
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namely if 


1 
2 

Ri .179 (60) 

Summing up, under this condition ^ is a contraction of Za,b,p into itself. Since ^(o)(?/) := 
2p^{y) = 2sin^ | belongs to Za^^^p, applying the Banach fixed point theorem we find 


7 < 




Theorem 2 Let 2;(o)(y) •=2sin^ Z(^n) ■=/^n :=with A^jl defined as in 

[43^ . The sequences {2:(n)}neN7 {/^n}neN converge respectively to the soliton solution z [in 
the norm !i44\ )] and to the corresponding for 7 at least in the range [60\) . With A 

defined as in (Ei), the errors of the n-th approximation are hound by 


A 




\hn-h\ < ^ 


32 L 


-77r 


-I A^ 


1-A 


11 ^( 1 )-^( 0)11 

(61) 


[To complete the proof we need just to note that, by (l57|) . the convergence of Z(^n) implies 
the convergence of fin and estimate the second error through standard arguments]. 

Remark 14]. 1 More refined computations of upper and lower bounds, with the present 
7-independent weight p^{y) = sin^ |, would show a 7-range of convergence of the above 
sequences slightly larger than (l60t) . By choosing a suitable 7-dependent weight p‘^{y), e.g. 
p'^iu) — ^(i)(y)/2) one could show that this range is actually significantly larger. This will 
be elaborated elsewhere. 


We explicitly work out the first approximation. We find: 


Zii){y) = v/l-7^2sin2 |-F7 


TT I COS ) +y — siny 


1 

hi = ^3 


(62) 

(63) 


6(1) (y) = 771- 



-|-const 


(64) 




1 

VT+TdaTvrTF 


7r7 

4a 


+ 0(7'). 


(65) 


The results are in good agreement with the numerical simulation plot in Fig. [2] (right). 
Note that the result (f65|) coincides with as announced. In a similar way one can 
determine iteratively solutions of type 3 (/r, z) even with low zm [he. not fulfilling the 
bound (|3] . 
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Appendix 

Classification of the singular points of llll|l . It is easy to check that the characteristic 
equations of m are 

+ = 0; (66) 

the upper, lower sign refer to any A^, respectively, 7 = 1 to Cfc. 

- The solutions Ai,A 2 for are both real, and distinct. is a saddle point and 
there are exactly fonr half-paths (called separatrices) with an endpoint on A^: the 
two ingoing represent motions approaching as ^ ^ 00 from the left or from the 
right, the two ontgoing represent motions leaving from as ^ — 00 towards the 
left or towards the right. 

- The solutions Ai, A 2 for are: 

— Both real if > 2 ( 1 — 7 ^)^/^. B^ is a node, and there are an infinite nnmber of 
half-paths ingoing to Bj. with the same tangent. These represent overdamped 
motions ending in B^ as ^ ^ 00 . 

— Complex conjngates (but not purely imaginary) if 0 < /r < 2(1 — 7 ^)^/^. Bk is 
a focus, and there are an infinite nnmber of half-paths ingoing to B^ along a 
spiral. These represent damped oscillatory motions ending in Bj. as ^ ^ 00 . 

— Opposite imaginary if fj. = 0. Any Bk is a center, and there exist closed paths 
(cycles) encircling it. These represent periodic motions around Bk, i.e. periodic 
oscillatory motion of the particle around g™. 

- If 7 = 1 Ai = 0, A 2 = —/r and if /i > 0 Cfc is a saddle-node: there only two half-paths 

(separatrices) in the half-plane g > {2k + l/ 2 )' 7 r (the ingoing represents a motion 
approaching as ^ ^ 00 from the right, the outgoing a motion leaving from Ak 
towards the right as ^ > — 00 ) and infinitely many in the half-plane g< {2k + l/ 2 ) 7 r 

(overdamped motions coming from the left and ending in Ck)- 

Proof of Prop. [H Let 0 < 20,2 <^ 0 , 1 , Zj{g) := z{g-, go, zqj] {j = 1,2) be the 
corresponding solutions of (jl3h and Gj the corresponding intervals giving their (maximal) 
domains. By continuity the inequality 


- -22 > 0 (67) 

will hold in a neighbourhood of go within Gi H G2. In fact, it will hold for all g £ Gi H G2. 
If per absurdum this were not the case, denote by 5 G Gi n G2 the least g > go (resp. 
largest g<go) where Z1 — Z2 vanishes: Zi{g) — Z2{g) = 0 ; then the problem (fTSll with initial 
(resp. final) condition z{g) = zi{g) = Z2{g) would admit the two different solntions zi,Z2, 
against the existence and uniqueness theorem. As for the monotonicity of g±, by the same 
theorem 2;i(5'2±) >-22(fi'2±) = 0 implies gi^>g2+ if <72-1-<00, otherwise <71+=52+= 00, and 
gi- < g2- if 52- > —00, otherwise gi- = 52- = —00. 

Proof of Prop. [2l Let 7^1 < /i 2 , 7ie > 72 ^, with one of the two inequalities being strict; 
for j = 1,2 let Uj{g) := u{g; go,uo', fJ-j,"/) be the corresponding solutions of (fT 3 |) with the 
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same condition Uj{go) = uq, and Gj the intervals giving their (maximal) domains, 
find 


72 - sin^f , 7i - sm^ 

U2g = -M2 H-7^- < + 


U2 


U2 


By the comparison principlqj (see e.g. [25]) it follows, as claimed, 


We 


ui{g)>U2{g) ge\gQ,g+[, ui{g)<U2{g) ge\g-,gQ[. ( 68 ) 

If e > 0, this implies: uiid) > hmgjg 2 + ^ 2 ( 5 ) = 0 and therefore 51 + > 52 + (the 

inequalities being strict as long as < 72 + < 00 ); limg|g^_ U 2 > 0 and therefore < 71 ^- > g 2 - (the 
inequalities being strict as long as gi- > — 00 ). Moreover, let gj{C) = 9 {^] fij,'yj) 

be the corresponding two solutions of (fH|l . i.e. the solutions of Q. We find 


52(0 = ^2(52(0) 


< ui{g2{0), 

> Uli92i0), 


Ve > ?o, 

ve < 6, 


while 52 (^ 0 ) = 9o = 5i(?o)- By the comparison principle this implies as claimed 92(0 < 
9 i{C) for all ^ € Xi n X 2 . Similarly one argues if e < 0. 

Proof of Prop. [3] Consider the Cauchy problem (|13p in subsequent intervals 
]gk,9k+i['^ G. Since the equation is invariant under g —> g + 2'K, by Prop. [T] if z{gi) 
is respectively larger, equal, smaller than z{go) then so are z{gk+i),Ik+i in comparison 
with z{gk),Ik respectively, for all k G K] in other words, the sequences {z{gk)}, {Ik} are 
either constant, or strictly monotonic. Eq. ([1^]) follows from (fT7)l applied in ]gk,9k+i[- 

If e = —, then rhs (|19l) > 27r7 > 0 for any k, so that the sequences are strictly increasing 
and diverging as A; —> 00 , whereas K must have a lower bound, otherwise z{gk) would 
become negative for sufficiently low k. 

If e = +, then the two terms at the rhs (|19p have opposite sign and can compensate 
each other. If the sequences are strictly increasing, the sides of ()19l) are positive for all k 
and Ik < 2'K’^lii. Applying (fTT]) to the interval [gk,gk+I^9\ for any Xg < 27r we find 

z{gk+I^ 9 ) - z{gk) = U{gk)-U{gk + Xg)-g j ds^2z{s). 

9k 

But \U{gk) — U{gk+Xg)\ is upper bounded, e.g. by 2+27r7, whence 


\z{gk+Xg) - z{gk)\ <2+27r7 + gh < 2 + 47 r 7 . (69) 

If per absurdum z{gk) diverged as /c —> 00 , then also z{gk+Xg) and in turn Ik [by (fTSl) ] 
would diverge, in contrast with Ik < 27r'y/g; so it must converge. Moreover, as before, K 
must have a lower bound. On the other hand, rewriting (II9h in the form z{gk-i) — z{gk) = 
glk-i — Svry, we see that if the sequences {z{gk)},{Ik} are strictly decreasing, the sides 
are positive for all k and larger than glQ — 2 tt^ > 0 for all negative k] this implies that 
they diverge as /c —> — 00 , and again by (l 6 ^ so do z{g), I{z, g). Whereas they must either 
converge as A: ^ 00 , or if must have an upper bound. 

"‘Here we recall the latter in the restricted version: if / fulfills conditions ensuring that the differential 
problem u' = f{x,u), u(xo) = u(xo), has a unique solution u, and u' < f{x,u) for all x, then it is u{x) < u{x) 
for all x>xo and u(x) > u(x) for all x<xo. 
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